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A coherent  optical  processor  for  displaying  a 
signal's  ambiguity  function  Is  described.  The  re- 
quired time  delay  Is  simulated  by  45  degree  rota- 
tions of  two  Identical  Input  transparencies,  and  the 
doppler  shift  by  a subsequent  one-dimensional  Fourier 
transformation.  The  entire  ambiguity  function  Is  dis- 
played In  the  output  (doppler  shift-time  delay)  plane. 
Examples  of  the  optically  computed  ambiguity  function 
for  single  and  double  pulse  signals  are  shown  to  be  In 
excellent  agreement  with  theory.  Advantages  of  this 
approach  over  other  schemes,  and  possible  extension  to 
real  time  processing,  are  also  discussed. 

Introduction 
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and  L2,  and  the  spatial  frequency  v Is  related  to  the 
horizontal  displacement,  x2,  on  plane  ?2  by9 

v - x2/xf  (3) 

By  an  appropriate  choice  of  an  Input,  this  coherent 
processor  will  be  shown  to  have  the  capability  of  ambi- 
guity function  display. 

Consider,  then,  the  one  dimensional  temporal 
signal  f(t)  In  Fig.  2a  and  Its  representation  on  the 
(t,v)  plane  [Fig. 2b],  By  rotating  this  function  coun- 
terclockwise about  the  origin  an  angle  of  e [F1g.2e], 
we  generate  the  function 


The  ambiguity  function,  first  Introduced  by  Wood-  f[tcose  ♦ tslne] 


i **ard  , has  been  applied  In  radar  In  predicting  the  ca- 
pability of  a given  signal  to  simultaneously  determine 
the  range  and  velocity  of  a target.  The  range  Is  de- 
termined by  the  time  delay,  t,  and  the  velocity  by  the 
doppler  shift,  v.  The  ambiguity  function  for  a given 
signal,  f(t),  is 

x(v,0  • £ f(t)f(t  - r)exp(-J2wvt)  dt  (1) 

In  optics,  Papoulls  has  employed  the  ambiguity  function 

In  analyzing  diffraction  phenomena  . 

In  this  paper,  we  describe  a rather  easily  Imple- 
mented coherent  processor  capable  of  generating  the 
afeblgulty  function  In  magnitude.  A similar,  yet  some- 
what more  elaborate,  scheme  for  generating  x(v,t)  In 
both  magnitude  and  phase  1$  also  presented. 

Such  a scheme.  Tor  example,  would  need  to  be  utilized 
when  further  coherent  processing  of  the  ambiguity  func- 
tion Is  required. 

Cutrona  et.al.2**  and  Preston®  have  proposed  a 

coherent  amblgutty  function  processor®  which  utilizes 
multiple  channels  to  display  the  ambiguity  function 
for  discrete  values  of  t.  The  scheme  of  Casasent  et. 

) al.2  generates  one-dimensional  ■slices''  of  the  ambigu- 
, Ity  function  In  the  (v,r)  plane.  Similar  one-dlmen- 
I slonal  displays  have  also  been  electronically 

^ produced8.  Our  method,  as  described  In  the  following 

sections,  (1)  displays  |x(v,t)|2  In  a continuous 
(rather  than  quantized)  form  over  the  entire  (v,t) 
plane,  (2)  has  the  capacity  for  extension  to  real  time 
processing,  and  (3)  Is  easily  Implemented. 

Implementation  Scheme 

The  coherent  processor  capable  of  displaying  the 
ambiguity  function  (In  magnitude)  Is  pictured  In  Fig  1. 
The  field  amplitude,  U(v,t),  in  plane  P{  Is  related 

to  the  coherently  Illuminated  transparency,  S(t,t)  In 
plane  Pj  by  a one-dimensional  Fourier  transform: 

U(v,t)  « exp(-J2*xf«2)  £ S(t,-r)exp(-J2fvt)  dt  (2) 

where  x Is  the  wavelength  of  the  spatially  coherent 
Illumination,  f Is  the  focal  length  of  both  lenses  L, 
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Thus,  for  a rotation  of  45  degrees,  we  obtain 
f((t  * t)/JZ),  and  for  a rotation  of  -45  degrees,  we 
obtain  f { (t  - r )/&).  If  we  cascade  transparencies  of 
these  two  functions.  In  the  Input  plane  of  the  proces- 
sor of  Fig. 1 so  as  to  form  the  product 
f[t  ♦ t)/v7]  f[t  - t)/v?],  then,  from  Eq.2,  the  resul- 
ting output  Is 

U(v,t)  • exp(-J2»xfv2) 

• C f |~~’j*xP(-J2»vt)  dt 

* exp[-J2wv(t  ♦ Xfv)]  (4) 

C T(t')f(t*  - /?r)exp[-J2»(^v)f]  dt' 

where  we  have  made  the  change  of  variable 


The  Intensity  distribution  associated  with  Eq.4  Is  Im- 
mediately recognized  as  a scaled  version  of  the 

squared  modulus  of  the  ambiguity  function:*8 
I(v,t)  - |U(v,t)|2 

• 2|X(^v./?t)  I*  (5 


Experimental  Results 

To  evaluate  the  performance  of  the  proposed  pro- 
cessor, the  ambiguity  functions  for  a single  and  doub- 
le pulse  signal  are  evaluated  analytically  and  com- 
pared to  the  corresponding  optical  systam  outputs.  In 
practice,  the  processor  output  Is  magnified  by  conven- 
tional means  for  observation  and  photographic  purposes. 

For  a Single  Pulse 

For  a single  pulse,  (F1q.3a),  we  may  write 

f(t)  - rect  (t/2T)  («) 

2 , ) Ibprprovnd  for  pmbllc  NlMN| 
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where  2T  Is  the  pulse  duration  and 


rect  (t) 


4 f 1 ; |t|  a 

l 0 ; t » 


The  geometric  Interpretation  of  f(t),  ft*  * *■),  and 

/ 7 

* T))f(tt  ~ Tl)  are  shown  in  Figs. 3b,  3c  and  3d 

/l  1 1 

respectively. 

Substituting  Eq.6  Into  Eq.1  followed  by  evaluation 
'yields  the  ambiguity  function 


The  corresponding  output  Intensity  Is 


|X(v.t)|  --  4(2T  - |t|)2Slnc2v(2T  - |t|) 

cos2(4wTv)  ; |t|  * 2T 
(2T  - |r|)2s1nc2v(2T  - |r|) 

; 2T  a | t | < 4T 
(6T  -|t|)  sine  v(6T  - |t|) 

; 4T  t |t|  ( 6T 
0 ; | t | t 6T 


. The  equations  descrlblnn  the  zero-value  loci  are  easily 

x(v,t)  -^(2T  - | t | ) s 1 nc  v(2T  - |t | )exp(-J.vT); |t  |*2T  thown  t0 

[o  ;|t|»2T  v ■ (2m  ♦ 1)/8T  ; |t|  t 2T 

<7)  v - n/(|t|  - 2T)  ; |t | < 4T  (13) 

where  sine  v 4 (sin  vv)/«v.  The  corresponding  output 

■ v • n/(6T  - » ) ; 4T  « |t  < 6T 

Intensity  Is 


v»t ) | * -^2T  - |t|)2$1nc2 


v(2T  - |t|)  ; | t | ( 2T 


For  purpose  of  Identification,  It  Is  Instructive 
to  examine  the  locus  of  points  where  the  ambiguity 
function  Is  Identically  zero.  From  Eq.8,th1s  zero  val- 
.ue  locus  may  easily  be  shown  to  be 


w ■ 7 r-  K '•  |t|  * 2T  (9) 

where  n Is  any  non-zero  Integer.  The  piecewise  hyper- 
bolic nature  of  these  curves  Is  shown  In  Fig. 4. 

The  ambiguity  function  for  a single  pulse  Is  gen- 
erated by  appropriately  rotating  two  Identical  thin 
silts  In  plane  of  the  coherent  optical  processor  of 

Fig. 2.  The  result  Is  shown  In  Fig. 5.  As  can  be  seen, 
the  coherent  processor  output  compares  quite  nicely 
with  the  theoretical  result  In  Fig. 4.  A three-dimen- 
sional computer  graph  of  the  corresponding  ambiguity 
function  modulus  may  be  found  In  Fig.  6.6  of  Rlhac- 

zek." 

For  a Double  Pulse 

For  a double  pulse,  (Fig. 6a),  we  write 
f(t)  • rect((t  ♦ 2T)/2T)  ♦ rect((t*2T)/2T)  (10) 

where,  for  convenience,  the  pulse  separation,  2T,  has 
been  chosen  to  be  equal  to  each  of  the  pulse  widths. 

The  geometrical  Interpretation  of 

1 Tl)  i,  shown  In  Fig. 6b.  The  ambiguity 

rt  ti 

function  associated  with  the  double  pulse  Is 

x(v.t)  -f  2(2T  - IxDsInc  v(2T  - ]t|) 

cos(4«Tv)exp(-j»Tv)  s |t|  < 2T 
-<2T  - | t | )s  1 nc  v(2T  - |t  J)exp(-J»Tv) 

; 2T  f t f 4T 

(6T  - 1 t | )s  1 nc  v(6T  - M)exp(-J.'tJ)  (11) 
; 4T  s t * 6T 
.0  5 It  » 6T 


where  m Is  any  Integer  and,  as  before,  n Is  any  non- 
zero Integer.  An  Illustration  of  these  zero-value  loci 
Is  presented  In  Fig. 7. 

By  appropriately  placing  two  Identical  double 
silts  In  planes  Pj  and  Pg  of  the  coherent  processor, 

the  ambiguity  function  for  the  double  pulse  Is  gener- 
ated. The  result,  shown  In  Fig. 8,  again  compares  quite 
favorably  with  the  theory. 

An  Alternate  Scheme 

The  coherent  processor  Just  discussed.  Is  capable 
of  displaying  the  ambiguity  function  only  In  magnitude. 
That  Is,  quadratic  and  linear  phase  factors  are  present 
on  the  output  plane  [see  Eq.4]. 

A processor  capable  of  generating  the  ambiguity 
function  both  In  magnitude  and  phase  is  presented  In 
Fig. 9.  In  plane  Pj,  we  place  the  transparency 

f{(T  ' t))  wh1ch  Is  formed  by  the  previously  discussed 

& 

45  degree  rotation  of  f(t),  followed  by  coordinate  re- 
versal (rotating  f((i.^..ll)  180  degrees  about  the  t and 

/z 

t axes).  The  scaling  lenses  and  have  respective 
focal  lengths  related  by 

f,  - ^ f2  (14) 

The  field  amplitude  Incident  on  the  left  of  plane  Pj  Is 

the  desired  f(t  - t).  This  will  multiply  the  transmit- 
tance f(t)  In  plane  Pj  to  give  Immediately  to  the  right 

of  P2  the  field  amplitude  f(t)f(t  - t).  The  geometrical 

interpretation  of  this  product  In  the  (t,0  plane,  for 
the  case  of  a single  pulse.  Is  shown  In  Fig. 10. 

With  reference  to  Eq.l,  It  remains  to  perform  a 
Fourier  transformation  with  respect  to  t.  This  Is  ac- 
complished with  cylindrical  lenses  La,  Lb,  and  lc  which 

have  respective  focal  lengths  of 

2f,  • fb  ■ 2fe  (15) 

One  sees  that,  from  plane  Pj  to  Pj,  Imaging  Is  perform- 
ed In  the  vertical  direction  by  la  and  lc  while  Fourier 


transformation  Is  Independently  performed  In  the  hori- 
zontal direction  by  L..  Thus,  the  field  amplitude, 
U(v,t),  1"  plane  Pj  Is  a scaled  version  of  the  ambigu- 
ity function: 

U(v,t)  ■ £ f(t)f(t  - T)exp(-J2*vt)  dt  (16) 

■ x(v,t) 

where  v Is  related  to  the  horizontal  displacement,  Xj, 
In  Pj  by 

v ■ x3/xf  (17) 

This  ambiguity  function  processor,  which  Is  more 
elaborate  than  that  shown  In  Fig. 2,  need  be  used  only 
when  the  ambiguity  function's  magnitude  and  phase  are 
required. 

As  a final  note,  we  observe  that  the  scheme  we 
propose  assumes  that  f(t)  Is  real  valued.  If,  In 
fact,  f(t)  Is  complex,  then  we  need  to  generate 
f*(t  - t) In  the  Integrand  of  Eq.l  In  order  to  properly 
generate  x(v,t).  Some  additional  work  appears  needed 
for  the  case  of  a general  complex-valued  f(t). 
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Fig.  1 . A coherent  processor  for  ambiguity  function 

display.  Both  the  lenses  have  focal  length  f. 
Fourier  transformation  Is  performed  In  the 
horizontal  direction  and  Imaging  In  the  verti- 
cal direction. 
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c1g.2.  (a)  A function  f(t)  In  time  and  (b)  on  the 

(t,r)  plane,  (c)  By  rotating  f(t)  counterclock- 
wise an  angle  of  a about  the  origin  of  the 
(t,t)  plane,  we  generate  f(tcese  ♦ ulna),  (d) 

The  function  f |n  the  (t.t) 
plane.  1 1 
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Fig. 3.  A single  pulse,  (a)  In  time,  (b)  on  the  (t,») 
plane,  (c)  rotated  4b  degrees  on  the  (t,r) 

plane  to  form  ff^^  Tl/,  (d)  the  product  of 
two  pulses  rotated  45  degrees  and  -45  degrees 
on  the  (t,t)  plane  to  form  f(^£^  *)>. 
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Fig. 4.  Zero  locus  plot  of  the  ambiguity  function  of  a 
single  pulse. 


Fig. 5.  The  ambiguity  function  (modulus  squared)  dis- 
play for  a single  pulse,  as  generated  by  the 
coherent  processor  of  F1g.2. 


' "T  T ' 3T 

la) 


"'OS* 


Fig. 6.  (a)  A double  pulse,  (b)  The  corresponding  func- 
tion T^)f(^y  T)>  on  the  (t.t)  plane. 


Fig. 7.  Zero  locus  plot  of  the  ambiguity  function  of  a 
double  pulse. 


Fig. 8.  The  amolgulty  function  (modulus  squared)  dis- 
play for  a double  pulse  as  generated  by  the  co- 
herent processor  of  Fig. 2. 
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Fig. 9.  Coherent  processor  for  generating  the  ambigu- 
ity function  In  phase  and  amplitude.  Scaling 
lenses  l.  and  L.  have  a focal  length  relation 

fi  - . *v 


Fig. 10.  The  product  f(t)f(t  - t)  on  the  (t,t)  plane 
for  the  single  pulse  of  Fig. 9a. 
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